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Abstract 



We analyze the data of low-energy cosmic-ray p spectrum, recently pub- 
lished by the BESS Collaboration, in terms of newly calculated fluxes for 
secondary antiprotons and for a possible contribution of an exotic signal due 
to neutralino annihilation in the galactic halo. We single out the relevant su- 
persymmetric configurations and discuss their explorability with experiments 
of direct search for particle dark matter and at accelerators. We discuss how 
future measurements with the Alpha Magnetic Spectrometer (AMS) on the 
Shuttle flight may disentangle the possible neutralino-induced contribution 
from the secondary one. 
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I. INTRODUCTION 



A recent analysis M of the data collected by the balloon-borne BESS spectrometer on 
cosmic-ray antiprotons during its flight in 1995 (hereafter referred to as BESS95 data) has 
provided the most detailed information on the low-energy cosmic-ray p's spectrum currently 
available: 43 antiprotons have been detected, grouped in 5 narrow energy windows over the 
total kinetic-energy range 180 MeV < Tp < 1.4 GeV. With this experiment the total number 
of measured cosmic-ray antiprotons in balloon-borne detectors over a period of more than 
20 years has more than doubled. Most remarkably, the BESS95 data provide a very 
useful information over the low-energy part of the p flux, where a possible distortion of the 
spectrum expected for secondary p's (i.e., antiprotons created by interactions of primary 
cosmic-ray nuclei with the interstellar medium) may reveal the existence of cosmic-ray 
antiprotons of exotic origin (for instance, due to pair annihilation of relic particles in the 
galactic halo [j^-Q, to evaporation of primordial black holes PJTH1] or to cosmic strings [|TT|]). 
In fact, a possible discrimination between primary (exotic) and secondary p's is based on 
the different features of their low-energy spectra: in this energy regime (Tp < 1 GeV) 
interstellar (IS) secondary p spectrum is expected to drop off very markedly because of 
kinematical reasons |12j| , while exotic antiprotons show a milder fall off. However, as will 



be discussed later on, this discrimination power is somewhat hindered by solar modulation 
and by some other effects affecting particle diffusion in the Galaxy. 

In Fig. 1 we report the cosmic-ray p flux at the top of the atmosphere (hereafter 
referred to as TOA flux) measured by BESS95 F° r experimental data referring to other 
measurements with much less statistics see Refs. []2|-|| . Also displayed in Fig. 1 are the 
minimal, median and maximal fluxes expected for secondary antiprotons at the time of the 
BESS95 data taking. These fluxes have been derived with a procedure which is described 
in detail in Sees. II-V. 

A comparison of the BESS95 data with the theoretically expected fluxes for secondary 
p's, as displayed in Fig. 1, leads to the following considerations: i) the experimental data 
are consistent with the theoretically expected secondary flux, within the experimental errors 
and the theoretical uncertainties; however, ii) the experimental flux seems to be suggestive 
of a flatter behaviour, as compared to the one expected for secondaries p's. Thus, natural 
questions arise, such as: a) how much room for exotic p's would there be in the BESS95 
data, for instance in case the secondary flux is approximately given by the median estimate 
of Fig. 1, b) how consistent with the current theoretical models would be the interpretation 
of the BESS95 data in terms of a fractional presence of exotic antiprotons, and c) how this 
interpretation might be checked by means of independent experiments? In the present note 
we address these questions within an interpretation of a possible excess of p's at low energies 



in terms of primary antiprotons generated by relic neutralinos in the galactic halo [13 



The present analysis P4J is mostly meant to a clarification of many theoretical points 



which will be even more crucial, when a much more statistically significant experimental 
information on low-energy cosmic-ray antiprotons will be made available by forthcoming 
experiments: AMS on the precursor Shuttle flight in May 1998 and on the International 



Space Station Alpha (ISSA) in January 2002 ||15|| , the satellite-borne PAMELA experiment 
|T6|j and balloon-borne measurements |T7| . 

Our paper is organized as follows. In Sec. II we discuss the cosmic-ray IS proton 
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spectrum which will be subsequently employed in deriving the secondary antiprotons. In 
the same section we also illustrate how we treat the solar modulation to connect the IS 
spectra to the corresponding TOA fluxes. In Sec. Ill we discuss the sources of cosmic 
antiprotons, both of primary and of secondary origins. Cosmic rays diffusion properties are 
derived in Sec. IV; the TOA p spectra are given in Sec. V. In Sec. VI we compare our 
theoretical fluxes with the BESS95 data and single out the neutralino configurations which 
may be relevant for the present problem. Sees. VII and VIII are devoted to an analysis 
on how these supersymmetric configurations can be explored by direct searches for relic 
neutralinos and by experimental investigation at accelerators. Conclusions and perspectives 
in terms of the forthcoming measurements of low-energy cosmic-ray p's are illustrated in 
Sec. IX. 



II. COSMIC RAY PROTON SPECTRUM 

We have first to fix the primary IS cosmic-ray proton spectrum, since we need it for the 
evaluation of the secondary p's. The IS cosmic-ray proton spectrum is derived by assum- 
ing for it appropriate parametrizations and by fitting their corresponding solar-modulated 
expressions to the TOA experimental fluxes. 

Measurements of the TOA spectra have always suffered from large uncertainties, as 



discussed for instance in Ref. ||18|| . In the present paper we use the two most recent high- 
statistics measurements of the TOA proton spectrum: the one reported by the IMAX Col- 
laboration on the basis of a balloon flight in 1992 the other given by the CAPRICE 



Collaboration based on data collected during a balloon flight in 1994 |]20fl . These two fluxes 
are reported in Fig. 2. 

For the IS proton spectrum we have used two different parametrizations: one in terms 
of the total proton energy, E p = T p + m p , 

$J S (T P ) = Ap(E p /GeV)- - prQt ° nS (1) 
' nr • s • sr • GeV 

the other in terms of momentum, p (equivalent to rigidity for protons), 



=Bp-\p/GeVr 2 Pr0t ° nS r v , (2) 
v nr • s ■ sr ■ GeV 



where (3 = p/E p . For the solar modulation effect we have employed the Perko method |2T 
where the solar-modulated flux is given by 



The kinetic energies T and Tis are simply related by Tis = T + A, when T > T cut and 



by a more complicated relation otherwise [21 1. Thus, this solar-modulation recipe is fully 



defined, once the values for the two parameters A and T cut are given. 



The results of our best fits to the data of Refs. [|T^,pO] are reported in Table I in terms of 



the parameters of expressions ([!]) and (H) and of the solar-modulation parameter A. The first 
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and third sets of values for parameters A, a, A (for expression ([]])) and B, 7, A (for expression 
(g)) refer to 3-parameter fits over the entire energy range of the experimental data. These 
fits are mainly meant to fix the solar-modulation parameter A, since the low energy part of 
the spectra is strongly dependent on the effect of solar modulation. The second and fourth 
sets of values refer to 2-parameter fits (at fixed A) in the high energy range (T p > 20 GeV), 
where the solar modulation effect is less sizeable, however not negligible, and therefore the 
proper parameters of the IS flux (normalization and spectral index) can be determined more 
confidently. 

The best-fit values for T cut turn out to be always smaller than the value corresponding 
to the lowest T considered in the fit (i.e., T cut < 0.1 GeV). This is consistent with the 
determination of the cut-off rigidity of the diffusion coefficient in the heliosphere [f^,^3[ . 

From the values reported in Table I we notice that even using various parametric forms 
for the IS proton spectrum, the data of the two experiments of Refs. [|T9lpO| do not lead to 



a set of central values for the parameters mutually compatible within their uncertainties. 
This difference can be considered as due to systematics in the measurement of the proton 
spectra. We also notice that the parametrization of Eq.(Q) systematically provides larger 
values for the solar-modulation parameter A as compared to the ones obtained using the 
parametrization of Eq. ([!]). This is due to the steeper behaviour at low energies of the 
function of Eq. (0) with respect to Eq. (fj) . 

In Fig. 2 we display the curves of the best fits to the data of Refs. |l9| , |20l with the 



parametrization of Eq. (JJ). This is the expression for the IS proton spectrum that will be 
used hereafter. Thus, as a median cosmic-ray proton flux we take the expression of Eq.(p]), 
where the values assigned to the parameters A and a are the averages of the central values 
in the fits to the data of the two experiments for T p > 20 GeV, i.e. 

^f(T p ) = 15, 950 /3(iy GeV)-- 76 - prQt ° nS (4) 
' nr • s ■ sr • GeV 

To estimate the uncertainty in the cosmic-ray proton flux we have combined the uncer- 
tainties of the two parameters A and a (at fixed A) in our fits at high energies (second set 
of values in Table I). We found that a conservative (generous) uncertainty band is delim- 
ited by a minimal flux, given by expression ([I]) with values: A = 12, 300, a = 2.61, and a 
maximal flux given by expression (H) with values: A = 19, 600, a = 2.89. Minimal, median 
and maximal IS proton fluxes are displayed in Fig. 3 together with the experimental TOA 
spectra of Refs. . It turns out that, for instance, at 100 GeV the uncertainty in the 

IS proton flux, relative to the median spectrum, is < ±50%. 



III. PRODUCTION OF ANTIPROTONS IN THE GALAXY 



A. Secondaries p's 



Cosmic ray protons interact with the interstellar material that mostly spreads in the 
galactic disk. This conventional spallation is actually a background to an hypothetical 
super symmetric antiproton signal. It needs therefore to be carefully estimated, especially 
at low energies where new measurements are expected. The corresponding source term is 
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given by the convolution between the antiproton production cross section and the interstellar 
proton energy spectrum as 

5?*(r) = J E0 {^p ^ n n E p ) dE p , (5) 

where hh is the hydrogen density in the disk, v p the proton velocity and ipp is the proton 
density per energy bin at distance r from the galactic centre in the galactic frame. The 
collision takes place between an incoming high energy proton with an hydrogen atom at 
rest, lying in the gaseous HI and HII clouds of the galactic ridge. The proton energy is 
denoted by E p . It is larger than the threshold E® = 7 m. That spallation reaction may 
generate an antiproton with energy Ep. The relevant differential production cross section 
is the sum over the angle 9 between the incoming proton and the produced antiproton 
momentum 



— { E r E p } = 2, P t [ E,^ 



d{-cos9) , (6) 



LI 



where Pp = y Ep 2 — m 2 . That integral is carried out in the galactic frame, at fixed antiproton 
energy E p . The proton energy determines the center of mass frame (CMF) energy a/s = 



{2m (E p + m)} . The latter sets in turn the maximal CMF energy Ep mSLX which the 
antiproton may carry away once it is produced 



lm 2 



^pmax 2 ^ • ( 7 ) 

The range of angles 9 over which the integral in Eq. (|]) is performed is set by the requirement 
that the CMF antiproton energy Ep should not exceed the maximal value E^ max implied by 
kinematics. The Lorentz invariant antiproton production cross section Ep d 3 a/d 3 Pp has been 



parametrized by Tan and Ng |24J] as a function of the transverse and longitudinal antiproton 
CMF momenta Pp T and Pp L . We refer the interested reader to this analysis. The transverse 
momentum in the CMF is equal to Pp T = Ppsm8 while the longitudinal momentum P? L 
obtains from the component Pp cos 9 after a Lorentz boost from the galactic frame to the 
CMF of the reaction. Note finally that the antiproton production integral in Eq. (|5]) should 
be a priori performed everywhere in the confining magnetic fields of the galactic disk. It 
actually involves the interstellar proton flux <3> p which depends on the location r. 



B. p , s from neutralino annihilation 

The differential rate per unit volume and unit time for the production of p's from x~X 
annihilation is defined as 

lT Sy (Tp) - ^ = < <W, > 9(T P ) 2 , (8) 

where < cr ann t> > denotes the average over the galactic velocity distribution function of 
neutralino pair annihilation cross section cr ann multiplied by the relative velocity v of the 
annihilating particles, m x is the neutralino mass; g(Tp) denotes the p differential spectrum 
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u ami Uj± p p UJ -p 



( F) 

where F indicates the x~X annihilation final states, E> xh is the branching ratio into quarks 
or gluons h in the channel F and dNp/dTp is the differential energy distribution of the 
antiprotons generated by hadronization of quarks and gluons. In Eq.(H), Px( r ' z ) ls ^ ne mass 
distribution function of neutralinos in the galactic halo. Here we consider the possibility 
that the halo is spheroidal and we parameterize p x (r, z) as a function of the radial distance r 
from the galactic center in the galactic plane and of the vertical distance z from the galactic 
plane 

where a is the core radius of the halo, r is the distance of the Sun from the galactic center 
and / is a parameter which describes the flattening of the halo. Here we take the values: 
a = 3.5 kpc, r Q = 8 kpc. For /, which in principle may be in the range 0.1 < / < 1, we use 
the two representative values / = 0.5, 1 |p5fl . The quantity p° denotes the local value of the 
neutralino matter density. We factorize it as p x = £pi, where pi is the total local dark matter 
density. Here £ is evaluated as £ = min(l, Q x h 2 / (Qh 2 ) m i n ), where (fi/i 2 ) m ; n is a minimal value 
compatible with observational data and with large-scale structure calculations |2(J . All the 
results of this paper refer to the choice (flh 2 ) min = 0.03. The neutralino relic density fi x h 2 
is calculated as a function of the supersymmetric parameters as described in Ref. f27| . As 
for the value pi of the total dark matter density, this is calculated by taking into account 
the contribution given by the matter density of Eq. ([[(]) to the local rotational velocity. For 
instance, in the case of a spherical halo (/ = 1), a value of pi = 0.4 GeV cm -3 is obtained. 



When / < 1 (oblate spheroidal distribution), pi is given by p8| , p9 



M/)= W (/ = D /Ar f^L p . (ID 

All the quantities depending on the supersymmetric parameters have been calculated in 
the framework of the Minimal Supersymmetric extension of the Standard Model (MSSM) 
pOl , where the neutralino is defined as the lowest-mass linear superposition of photino (7), 
zino (Z) and the two higgsino states (H°, H%) 

X = 017 + a 2 Z + a 3 H° + aji^. (12) 

For the evaluation of the averaged annihilation cross section < a ann v > we have followed 
the procedure outlined in Ref. f|. We have considered all the tree-level diagrams which 
are responsible of neutralino annihilation and which are relevant to p production, namely: 
annihilation into quark-ant iquark pairs, into gauge bosons, into a Higgs boson pair and 
into a Higgs and a gauge boson. For each final state we have considered all the relevant 
Feynman diagrams, which involve the exchange of Higgs and Z bosons in the s-channel 
and the exchange of squarks, neutralinos and charginos in the t- and u-channels. Finally, 
we have included the one-loop diagrams which produce a two-gluon final state. For this 



annihilation channel, we have used the recent results of Ref. 31 
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The p differential distribution g(Tp) has been evaluated as discussed in Ref. ||. Here we 
only recall that we have calculated the branching ratios B xh for all annihilation final states 
which may produce p's, dividing these states into two categories: i) direct production of 
quarks and gluons, ii) generation of quarks through intermediate production of Higgs bosons, 
gauge bosons and t quark. In order to obtain the distributions dNp/dTp the hadronization 
of quarks and gluons has been evaluated by using the Monte Carlo code Jetset 7.2 [ 32| . For 
the top quark, we have considered it to decay before hadronization. 

We summarize now the main features of the MSSM scheme we employ here. The MSSM 
is defined at the electroweak scale as a straightforward supersymmetric extension of the 
Standard Model. The Higgs sector consists of two Higgs doublets Hi and H 2 , which define 
two free parameters: the ratio of the two vacuum expectation values tan f3 = (H 2 ) / (H\) and 
the mass of one of the three neutral physical Higgs fields; we choose as a free parameter 
the mass of the neutral pseudoscalar Higgs. The other parameters of the model are 
contained in the superpotential, which includes all the Yukawa couplings and the Higgs- 
mixing term /iHiH 2 , and in the soft-breaking Lagrangian, which includes the trilinear and 
bilinear breaking parameters and the soft gaugino and scalar mass terms. In order to deal 
with a manageable model, we impose the following usual relations among the parameters 
at the electroweak scale: i) all trilinear parameters are set to zero except those of the 
third family, which are unified to a common value A; ii) all squarks and sleptons soft-mass 
parameters are taken as degenerate: m z ~ = m$. = mo; iii) the gaugino masses are assumed 
to unify at Mqut, and this implies that the U(l) and SU(2) gaugino masses are related at 
the electroweak scale by M x = (5/3) tan 2 9wM 2 . 

When all these conditions are imposed, the supersymmetric parameter space 
is completely described by six independent parameters, which we choose to be: 
M 2 , /J*, tan/3, rriA, mo, A. In our analyses, we vary them in the following ranges: 10 GeV < 
M 2 < 500 GeV (21 steps over a linear grid); 10 GeV < < 500 GeV (21 steps, linear grid); 
75 GeV <m A < 500 GeV (15 steps, logarithmic grid); 100 GeV < m < 500 GeV (5 steps, 
linear grid); —3 < A < +3 (5 steps, linear grid); 1.01 < tan/3 < 50 (15 steps, logarithmic 
grid) . 

The supersymmetric parameter space is constrained by all the experimental limits ob- 
tained from accelerators on supersymmetric and Higgs searches. The latest LEP2 data on 
Higgs, neutralino, chargino and sfermion masses |33[] and the constraints due to the b — > s+7 
process |34 are imposed. Moreover, the request for the neutralino to be the Lightest Su- 
persymmetric Particle (LSP) implies that regions where the gluino or squarks or sleptons 
are lighter than the neutralino are excluded. A further constraint is imposed by requiring 
that all the supersymmetric configurations which provide a neutralino relic abundance are 
in accordance with the cosmological bound Q x h 2 < 0.7. 



IV. DIFFUSION OF COSMIC RAYS INSIDE THE GALAXY 

The propagation of cosmic rays inside the Galaxy has been considered in the framework 
of a two-zone diffusion model. We have followed here the same analysis as Webber, Lee and 
Gupta [p5|. The Milky Way is pictured as a thin disk, 200 pc across, that extends radially 



up to R = 20 kpc from the galactic center. That ridge lies between two extended layers 
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~ 3 kpc thick, where cosmic rays diffuse in erratic magnetic fields. Mere diffusion governs 
the propagation of the particles in the disk and in the confinement regions that extend on 
either side. Assuming that steady state holds, the proton density i/j p , per energy bin, at 
some location r and z, is given by 

^ = = V-(iTV^p) + 2h5(z)q(r) - 2h5(z)T p ^ p . (13) 

The diffusion coefficient K is assumed to be essentially independent of the nature of the 
species that propagate throughout the Galaxy. It increases with rigidity 1Z according to the 
relation 



K 



0.6 



Kill) = K (l + -J , (14) 

where Kq = 6 x 10 27 cm 2 s -1 and TZq = 1 GV. Below that critical value, the diffusion 
coefficient stays constant while above 1 GV, it increases like TZ 06 . Sources are located 
in the galactic ridge at z — 0. Their radial profile is inferred from the survey by Lyne, 
Manchester and Taylor of the galactic distribution of stellar remnants and pulsars [^6j with 
q(r, 0) oc p a exp(-bp) where p = r/R, a = 0.6 and 6 = 3. Finally, cosmic ray protons may 
interact with the interstellar gas. The latter is assumed to be concentrated in the disk. The 
probability per unit time that a proton collides with an interstellar hydrogen atom at rest is 

T p = n H agt/p • (15) 

The hydrogen density Tin is assumed to be constant all over the disk. The value of nn = 1 
cm -3 is basically consistent with measurements of the hydrogen column density derived from 
HI and CO surveys. It implies in particular a maximal value of ~ 9 x 10 22 H cm -2 to be 
compared to an average of 5 x 10 22 H cm -2 on the observations of the galactic center. The 
densest spot is inferred from CO measurements to reach a level of ~ 1.4 x 10 23 H cm -2 . 
The total interaction cross section between the propagating high energy protons and 
the hydrogen atoms of the interstellar medium has been borrowed from the work by Tan 
and Ng [j3^] . Above a kinetic energy of 3 GeV, it may be expressed as 



= (32.2 mb) {1 + 0.0273 U} , (16) 

where the parameter U is defined as 

U = In (£ p /200 GeV) . (17) 

Below T p = 3 GeV, expression ( |16"D needs to be divided by a low energy correction factor 
equal to 1 + 0.00262 T p ~ Cp where 

C p = 17.9 + 13.8 In T p + 4.41 In 2 T p . (18) 

The galactic disk is assumed to be infinitely thin, hence the factor 2h5(z) in the diffusion 
equation flT3|), where 2h = 200 pc stands for the actual thickness of the ridge. 

Because the Galaxy is axi-symmetric, we can expand the proton density ip p as a series 
of Bessel functions of zeroth order 
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+oo 
i=l 



(19) 



where p = r/R, while on is the i-th zero of the Bessel function Jo. The proton density is 
ensured to vanish at the radial boundary r = R of the system. The Bessel transforms Vi 
must also drop to zero at the boundaries of the confinement regions, at a distance L = 3 
kpc from either side of the galactic disk. The distribution of cosmic ray sources may also be 
expanded as a series of Bessel functions 

+ °° a- 

q(r,E p ) = J2frQtot(E p )J (a iP ) , (20) 
i=i zn 

where Qtat(E p ) stands for the total galactic rate of production, per energy bin, of cosmic 
ray protons with energy E p . The Bessel transforms qi are readily inferred from the radial 
distribution of the sources in the galactic disk 

"< = ^ jfcT) U ' q{p) dp 1 {/»' q(p) dp2 } " ' (21) 

Bessel expanding the diffusion equation flT3] ) leads to simple differential relations which 
the functions Vi(z) satisfy. The latter are even functions of the height z that vanish at the 
boundaries of the diffusion layers. A straightforward algebra leads to 

nz,E p ) = ^Q tot (E p )sinh{|(L-|z|)}/sinh{^} , (22) 

where Si = 2ai/R and where the coefficients A4 are defined by 

Ai = 2hT p + KSiCoth^Y^j . (23) 

Because the diffusion term dominates the behaviour of the coefficients Ai, the proton 
energy spectrum does not vary much all over the Galaxy, except for a global normalization 
factor. In other words, the ratio of the proton fluxes taken at two different energies is quite 
insensitive to the location M, hence 

® p (M,E 1 )/<S> p (M,E 2 ) ~ V i (0,E 1 )/V i (0,E 2 ) . (24) 

This will turn out to be important when we compute the energy spectrum of secondary 
antiprotons. 

The two-zone model is a refinement with respect to the old leaky box scheme. The 
confinement layers are necessary in order to account for the low abundance of the 10 Be 
unstable element with respect to its stable partner 9 Be. The former nucleus has a half- 
lifetime of 1.6 million years (My) and plays the role of a chronometer. Observations indicate 
that cosmic rays are trapped in the magnetic fields of our Galaxy for approximately 100 
My before they escape in the intergalactic medium. On the other hand, the amount of 
secondary light nuclei such as lithium, beryllium and boron (Li-Be-B) is well explained by 
the spallation of primary carbon, oxygen and nitrogen (CNO) nuclei. The latter spend a 
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mere 5 My in the galactic plane where they cross a column density of ~ 10 gcm~ 2 . Cosmic 
rays are therefore confined most of the time in extended reservoirs above and beneath the 
matter ridge, where they just diffuse without interacting much with the scarce interstellar 
medium. We have estimated the grammage which the CNO elements cross during their 
journey inside the galactic disk. Their distribution is inferred in just the same way as for 
the protons. The average electric charge per nucleon is now 1/2 instead of 1 for the protons, 
hence a slightly modified relationship between the kinetic energy per nucleon and the rigidity 
of the nucleus under consideration. The grammage is defined as the product 

A e = n H v N r disk , (25) 

where the confinement time in the disk alone is denoted by r^k- The escape length A e is 
expressed in units of g cm -2 . Because cosmic rays either escape from the disk or interact 
with its gas, the total number J\f$ lsk of particles contained in the galactic ridge satisfies the 
balance relation 

A /"disk 

Qn = ^- + r N < isk . (26) 

Tdisk 

The rate Qn at which the CNO primaries are produced is set equal to the sum of the escape 
rate from the galactic ridge and of the interaction rate with the interstellar gas. Notice that 
in the case of the two-zone model, the amount A/^ lsk of cosmic rays travelling in the disk 
alone may be expressed as the series 

^=4./^g ^AM. (27) 

In the coefficients A^, the relevant cross section that accounts for the interactions of the 
CNO species with the interstellar hydrogen has been averaged at a mere 250 mb. In Fig. 
4, an estimate of the grammage A e crossed by the CNO elements is presented as a function 
of the kinetic energy per nucleon (solid line). It reaches a maximum of ~ 8 g cm -2 at 500 
MeV/n. It decreases at low energies with the velocity. It also drops at high energies as a 
result of a better diffusion and hence a lower time of residence in the disk. The dashed curve 
refers to the grammage of the protons. At fixed kinetic energy, the diffusion coefficient is 
slightly smaller for these species than for heavier elements, hence a larger escape length A e . 
Measurements of the 2 H abundance have been performed |38j from the Voyager probe at a 
distance of 23 AU and at energies lying between 20 and 50 MeV/n. With a solar modulation 
parameter of ~ 360 MV, this translates into an energy of ~ 230 MeV/n in interstellar space. 
The analysis by Seo et al. of these data is well accounted for by the leaky box model using 
a grammage A e (B/C) ~ 8 g cm" 2 . This is in excellent agreement with the results of our 
two-zone model presented in Fig. 4, where the diffusion coefficient K is given by relation 
(|i~4]). Ficenec et al. |39j have taken data on 3 He between 100 MeV/n and 1.6 GeV/n. They 



conclude that the grammage of primary cosmic rays is well fitted by A e = (10.5 + 2.5 — 2.8) f3 
g cm -2 . The extreme values of that fit are featured by the dotted curves of Fig. 4. Notice 
that the CNO grammage inferred from our two-zone model lies in the range of escape length 
delineated by the Ficenec et al. extreme values, for energies in interstellar space comprised 
between 200 MeV/n and 1.5 GeV/n. The expression which we have adopted for the diffusion 



10 



coefficient K is therefore well supported by measurements j|(J of the grammage encountered 



by primary CNO cosmic rays while they propagate within the galactic ridge. 

The propagation of antiprotons throughout the Galaxy follows the same trends as for 
the protons. We focus first on the species produced by the spallation of cosmic ray protons 
with the interstellar gas of the ridge. Their density tpp, per energy bin, follows the diffusion 
equation 

V-(^V^) - 2h5(z)T p ^ p + 2h5{z)qf\r) - 2h5(z) A {b{E)%} = , (28) 

where steady state has once again been assumed. We recognize the usual diffusion term 
as well as the contribution due to the interactions of the antiproton cosmic rays with the 
matter of the disk. The total interaction cross section between antiprotons and the hydrogen 
atoms of the interstellar medium has also been borrowed from the analysis by Tan and Ng 
P7| . Above Tp = 50 MeV, it may be parametrized as 

<7*g = (24.7 mb) {l + 0.584 TjT ai15 + 0.856 Tr - 566 } , (29) 

where the antiproton kinetic energy Tp is expressed in GeV. The spallation source term has 
already been discussed in Sec. Ill A. It obtains from the convolution @ of the antipro- 
ton production cross section with the proton energy spectrum. In order to simplify the 
calculations, we define the effective antiproton multiplicity 



™ = « k {E * - E > } m dE ° ■ (3o) 

Because the ratio of the proton fluxes taken at two different energies does not depend on 
the location, the effective antiproton multiplicity N? s is inferred to be only sensitive to the 
energy. It is therefore constant throughout the galactic ridge and may be computed once 
and for all as a function of the energy Ep of the produced antiproton before the diffusion 
equation ( |2"8"| ) is solved. The spallation production term readily simplifies into 

qf* k (r, E) = <j$(E) Nf(E) v p n H ^{r, E) . (31) 

Under that form, it may be immediately expanded as the usual series of Bessel functions 
of zeroth order. The last term in relation fl28|) stands for the energy losses suffered by 
the antiproton cosmic rays while they propagate in the galactic disk. That term actually 
exists for any cosmic ray species. Because the particle fluxes do not significantly drop at low 
energies, this effect is in general neglected. Fluxes tend even to increase below 1 GeV. In the 
specific case of secondary antiprotons, that is no longer valid. Because a high energy proton 
has very little chance to produce an antiproton at rest while colliding on an hydrogen atom, 
the secondary antiproton flux sharply drops when the energy decreases below ~ 1 GeV. 
Energy losses tend to shift the antiproton spectrum towards lower energies with the effect of 
replenishing the low energy tail with the more abundant species which had initially a higher 
energy. This process is understood here as a mere diffusion in energy space. The rate at 
which the antiproton energy varies b(Ep) = Ep takes into account two main effects. First, 
antiprotons may suffer from ionization losses while they travel across the interstellar gas. 
This mechanism yields the following contribution to the energy loss rate 
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bUE) = -iirrl m e c 2 n H °- {in (^~) + In (>V) - P'j ■ (32) 

In molecular hydrogen, the ionization energy E has been set equal to 19.2 eV; here 7 = E/m. 
The classical radius of the electron is denoted by r e and the electron mass is m e . Second, 
the dominant contribution to the energy losses arises from the elastic scatterings of high 
energy antiprotons on the hydrogen atoms of the disk. This mechanism is a counterpart to 
the collision process whose rate is Tp. An antiproton with initial energy E\ ends up after 
such a collision in a final state with the lesser energy E 2 . Elastic scatterings feed therefore 
the low energy part of the antiproton distribution. They have been described here as if 
they induced a continuous change in the antiproton energy. Our assumption is correct on 
average, hence the contribution 

&scat(£p) = - y {af u (E p )n u v p } . (33) 

The elastic cross section cr?j^ obtains from the difference crp^ — where the annihilation 
cross section is given by 

a™ = (661 mb) {l + 0.0115 Tr°- 7U - 0.948 T p a0151 } , (34) 

between 100 MeV and 12 GeV, i.e., the energy range under scrutiny here. Low energy data 
are fairly consistent with an average energy loss approximately equal to a half of the initial 
antiproton kinetic energy . 



The antiproton density ipp, per energy bin, may be Bessel transformed into the functions 
Vi whose variations with height z are given by 

Pi(E»z) = ^(^p,0)sinh{|(L-|z|)}/sinh{^} . (35) 

In the galactic disk at z = 0, the Bessel transforms Vi(Ep, 0) only depend on the antiproton 
energy E p . They actually satisfy a first order differential equation 

2h^{bV l ) + Bfc = 2h (o$Nfv p )\ E n H Vi(E,0) , (36) 

which we have numerically solved for each order i < 100. At high energy, antiprotons are 
insensitive to the energy losses. Starting therefore from an unperturbed spectrum, we have 
decreased the kinetic energy from 10 GeV down to 100 MeV while integrating equation (BBI). 
The coefficients Bi obtain from A4 by replacing the rate T p by its antiproton counterpart Tp. 
The above mentioned method has been applied to the case of the median proton flux (f|) 
derived from the IMAX and CAPRICE measurements. The solid curve of Fig. 5 stands for 
the corresponding antiproton interstellar flux. Energy losses have been taken into account. 
This is not the case however for the dot-dashed line where the same proton spectrum has 
been assumed. Note that energy losses tend actually to replenish the low energy part of 
the antiproton distribution. This effect is particularly evident at low energy. For Tp ~ 100 
MeV, the antiproton flux increases by more than an order of magnitude when energy losses 
in the gaseous disk are considered. At larger energies, the upward shift of the spectrum is 
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less sizeable. For an interstellar kinetic energy of 600 MeV, the increase has reduced to ~ 
30%. Even in the case of minimal solar modulation, 600 MeV in interstellar space translate 
into Tp ~ 300 MeV. The dotted and dashed curves respectively stand for the antiproton 
spectra derived from the minimal and maximal IS proton flux discussed at the end of Sec. II. 

We finally discuss the case of the antiprotons produced in the annihilation of neutralinos 
potentially concealed in the galactic halo. The diffusion equation is quite similar to relation 



V-(^V^ P ) - 2h8(z)T p iP P + qr y (r,z) = . (37) 

Because the energy distribution of these supersymmetric antiprotons is fairly flat, energy 
losses in the disk should play a negligible role. They have not been considered here. The 
source term (|8]) has already been discussed in Sec. Ill B. The antiproton production extends 
now all over the Galaxy and not solely in the disk. The solution of the diffusion equation 
(|37j) follows however the same trends as for the previous cases. The antiproton energy 
distribution ipp may still be expanded as a series of its Bessel transforms Vi(Ep,z). Since 
energy losses are negligible, the latter obey the simple differential equation 

K " - 2M(z)r p ^ + C Sy {z) = . (38) 



The Bessel transforms of the supersymmetric antiproton source distribution q S p hy are defined 



as 



1 f 1 

<?r Sy W = nf \ / J o( a iP)% Sy ( r = pR,z)dp 2 . (39) 
Ji{a i ) Jo 



Outside the galactic ridge, equation fl38f) simplifies even further into 

d 2 V- a 2 a susy (z) 

111 _ 2b. p. + 1l \21 = n (40) 

dz* R 2 K U ' 1 Uj 

The general solution may be expressed as 

Pi(z) = a cosh + h sinh (^f) (41) 



2 J V 2 



L exp(-^|z-z'|/2) qr y (z') dz' 



KSi Jo 

where Si = 2a.il R. We leave as an exercise the determination of the constants of integration 
a,i and They obtain from the requirement that the Bessel transforms Pi vanish at the 
boundaries z = ±L of the confinements regions that extend on either side of the ridge. 
Because the antiproton distribution is an even function of the height z, the interested reader 
may also show that Vi(0) = hTpVi(0) / K . The final result readily obtains as 

*w = i^{*w ■ < 42 > 
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This expression describes the actual propagation of antiprotons which have been produced 
in remote regions of the halo and that propagate backwards in the magnetic fields of the 
Galaxy. The functions T{z) and Q{z) are respectively defined by 

T{z) = [ Z sinh ( ^\z - z'\) qr y (z')dz' , (43) 

and 



o 



G(z) = 2hT p sinh (^) + KS t cosh . (44) 

The interstellar flux at the solar system of the antiprotons produced by the annihilation 
of hypothetical supersymmetric species comprising part of the galactic halo may now be 
expressed as 

$p(©,T p ) = <a aan v> g(T p ) C susy (T p J). (45) 

I m x J 

The density of reference po has been set equal to 1 GeV cm -3 . The coefficient C susy (Tp, f) 
is defined as 

C susy (T p J) = i-^ P efr (©,Tp) . (46) 

The effective energy distribution ipp cS is taken at the solar circle and has been derived with 
the above mentioned method where an effective antiproton source term {p x (r, z) / po} 2 has 
been assumed. The latter term depends on the flattening / of the halo. Note that C susy 
is not a flux of particles. It is a mere coefficient that is actually expressed in units of cm 
sr _1 . Fig. 6 illustrates the behaviour of this coefficient when the antiproton kinetic energy 
is varied from 100 MeV up to 10 GeV, for three different values of the flattening factor / = 
0.1, 0.5 and 1. The coefficient C susy (Tp, /) exhibits a smooth maximum around T p ~ 1 GeV. 
Below that value, it tends to decrease with the antiproton velocity like Vp/Bi. For higher 
energies, the diffusion takes place more efficiently and the cosmic rays escape more easily 
from the galactic magnetic fields, hence a lower density in the disk. When the flattening 
increases, the dark matter halo is compressed towards the ridge. There are many more 
neutralinos in the diffusion layers where antiprotons are kept confined, hence a larger flux. 
The evaluation of the IS p flux due to neutralino annihilation is then performed by using 
Eq.(iU). 



V. PRIMARY AND SECONDARY ANTIPROTON TOA FLUXES 



Our TOA antiproton fluxes are derived from the corresponding IS spectra, by employing 
the Perko solar-modulation procedure | 21|j , already defined in Sec. II. In that section we also 
derived the values for the parameter A relevant to the measurements of Refs. [ |T5||Z0"| ] . In order 
to obtain the A values to be applied in case of experiments performed at different times, we 
use the results of Papini, Grimani and Stephens (PGS) [Q. These authors derived simple 
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analytic expressions as best fits to the measured spectra of the TOA primary cosmic-ray 
protons, obtained from a large collection of data over a couple of solar cycles. They provide 
the parameters of these fits for periods of maximum and minimum solar activity. By fitting 
their analytic expressions with the solar-modulated flux derived from our parametric form 
of Eq. ([[]), we find the following average values for A at minima and maxima: A m j n = 320 
MeV and A max = 800 MeV, respectively. 

In Fig. 7 we plot the time variation of the solar-modulation parameter A, as obtained 
by our best fit to the experimental data. The full circles represent the best-fit values to 
the PGS average fluxes at minima (MIN) and at maxima (MAX) and to the fluxes of Refs. 
|T9| , |20[| . The open circle refers to the BESS95 data taking period. The cross denotes the 
extrapolated value at the time relevant for the AMS measurements with the Shuttle flight. 

In Figs. 8 and 9 we display how the effects of the flux distorsion at low energies, induced 
by solar modulation, is much stronger for the primary flux than for the secondary one. This 
may be simply understood in terms of the nature of Eq.(^) and of the different shapes for 
secondaries and primaries. 



VI. COMPARISON WITH EXPERIMENTAL DATA 

Early measurements of cosmic-ray antiprotons have been plagued by low-statistics prob- 
lems and brought to serious conflicting results at low energies (Tp < 0.4 GeV) in the past 
43| . As already mentioned in our introduction, a recent analysis [|1| of the data collected by 



the BESS spectrometer during its 1995 flight (BESS95) has provided a significant improve- 
ment in statistics in the low-energy region, with a total of 43 p's in the kinetic-energy range 



180 MeV <Tp < 1.4 GeV ||44|| . This then allows an interpretation of the experimental data 
in terms of theoretical models in a more meaningful way than in the past. A further substan- 
tial breakthrough in this direction will be provided by the forthcoming measurements with 
AMS lfL5|l , the satellite-borne PAMELA experiment |16| and balloon-borne measurements 
IT! 



The BESS95 data are displayed in Fig. 1 and compared with our theoretical evaluations 
for secondary antiprotons. Our curves are derived according to the procedure outlined in 
previous sections. Solar modulation is evaluated at the time of the BESS95 measurement. 
The band delimited by the dotted and the dashed curves provides the uncertainty in the 
secondary p's flux due to the corresponding uncertainty in the primary IS cosmic ray proton 
flux (see Sec. II). It turns out that this uncertainty is < ±30% for Tp < 2 GeV and it grows 
up to ±50% at Tp ~ 10 GeV. 

From a first look at Fig. 1 it is apparent that the experimental data are consistent with 
the flux due to secondary p's. This is indeed quantitatively confirmed by a x 2 ~evaluation, 
which shows that our median curve for secondaries fits the BESS95 data with a low reduced- 
X 2 value: (x 2 ) red = 0.83 (for 5 d.o.f.) |45 . 



However, it is interesting to explore which would be the chances for a signal, due to relic 
neutralino annihilations, of showing up in the low-energy window (Tp < 1 GeV). This point 
is very challenging, especially in view of the interplay which might occur among low-energy 
measurements of cosmic-ray p's and other searches, of quite a different nature, for relic 
neutralinos in our Galaxy. 



15 



Actually, we find that the agreement between BESS95 experimental data and theory 
may be improved by adding a fraction of neutralino-induced p's to the standard secondary 
antiprotons. The best fit to the experimental data with a total theoretical flux $ th = $^ed + 
$ susy , performed by varying the supersymmetric parameters over the grid defined in Sec. IIB, 
provides a value (x 2 )red = 0.28, with an improvement over the (x 2 )red previously obtained 
by using the secondary flux only. This fact cannot certainly be taken as significant of an 
evidence of a neutralino-induced antiproton signal, but shows that indeed the low-energy 
region p spectrum is still a quite interesting window for exploring p's of supersymmetric 
origin, and encourages further investigation of the problem. 

Now we wish to specifically determine which regions of the supersymmetric parameter 
space (and then which neutralino configurations) might be relevant for the problem at hand 
and how these could be investigated by other experimental means. As a quantitative cri- 
terion to select the relevant supersymmetric configurations, we choose to pick up only the 
configurations which meet the following requirements: i) they generate a total theoretical 
flux $ th which is at least at the level of the experimental value (within 1-a) in the first 
energy bin; ii) their (x 2 )rcd, in the best fit of the BESS95 data, is bounded by (x 2 )rcd < 
2.2 (corresponding to 95% C. L. for 5 d.o.f.). This set of configurations is hereafter de- 
noted as set M; its subset, whose Q x h 2 values fall in the cosmologically interesting range 
0.03 < Q x h 2 < 0.7, is denoted as set N. An example of a fit to the BESS95 data which 
includes a neutralino-induced signal with a (x 2 ) rc d < 2.2 is shown in Fig. 10. This signal 
corresponds to a neutralino with the following properties: m x = 62 GeV, P = 0.98 and 
tt x h 2 = 0.ll. 

On the other hand, supersymmetric configurations with a (x 2 )rcd > 4 have to be consid- 
ered strongly disfavoured by BESS95 data (actually, they are excluded at 99.9 % C.L.). We 
call R this set of supersymmetric configurations and we will discuss them later on. Super- 
symmetric configurations belonging neither to M nor to R can only provide a p flux fully 
buried in the secondary p flux and are then completely irrelevant for the problem under 
discussion. 

The composition of configurations in sets M and N are displayed in Fig. 11 (Fig. 12), 
where tan j3 (m x ) is plotted in terms of the fractional amount of gaugino fields, P = al+a 2 ,, in 
the neutralino mass eigenstate. From Fig. 11 we notice that gaugino configurations are more 
numerous than others, with only a slight correlation with tan (3; the requirement of a sizeable 
contribution to the relic density introduces a noticeable reduction in the number of higgsino- 
like and mixed configurations. Fig. 12 shows that higgsino-like and mixed configurations 
are much stronger constrained in the neutralino mass range than the gaugino-like ones. In 
Fig. 13 we display the features of configurations of set R. These configurations, which are 
to be considered excluded on the basis of the BESS95 data, turn out to be gaugino-like with 
masses on the low side. 

Up to now, we have discussed our results in terms of a spherically symmetric galactic 
halo. The effect of a flattening in the dark matter distribution is to enhance the primary 
p flux. Since the size of this flux is proportional to the function C susy (Tp, f), defined in 
Sec. IV, the enhancement of the primary flux as a function of / may be read directly from 
Fig. 6. For instance, for / = 0.5 the enhancement factor is 2.3. This has consequences 
on the nature of configurations in sets M, N and R. By way of example, we plot in Fig. 
14 the scatter plot for configurations of set R for a flattening of / = 0.5. This may be 
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compared with the corresponding plots of Fig. 13 which refer to / = 1. Obviously, the 
enhancement of the primary flux, induced by the halo flattening, increases the number of 
excluded configurations. 

In the present paper we have considered only uniform dark matter distribution inside 
the density profile of Eq. (PX|). As is well known, any effect of local density enhancement or 
dumpiness would induce a substantial increase in the primary p spectrum, as in any other 
signal due to pair annihilation taking place in the halo [J5|. 

Let us now examine whether our relevant neutralino configurations may be explored in 
terms of direct detection experiments for particle dark matter candidates. 



VII. EXPLORATION BY DIRECT DETECTION OF RELIC PARTICLES 

The measurements of the energy differential rates in experiments of direct search for 
particle dark matter enable the extraction of an upper bound for the neutralino-nucleon 
scalar cross-section o"° c ^ r on , multiplied by the neutralino local (solar neighbourhood) density, 
i.e. an upper bound for the quantity £cr" c ^° n , once a specific value to the total local dark 



matter density is assigned ||47|| . By combining all present experimental data |4q] , we obtain 
the (90 % C.L.) upper bound displayed in Fig. 15 by the open solid curve (the total 
local dark matter density is normalized here and in the rest of this paper to the value 
pi = 0.4 GeV cm -3 ). The experiments which are essential in the determination of this 



upper bounds, in the neutralino mass range considered here, are those of Refs. j4"9" |. The 
region in Fig. 15 delimited by a closed contour is the one singled out by the experiment 
of Ref. |)U| as possibly indicative of an annual modulation effect (for an interpretation 
of these data in terms of relic neutralinos see Ref. |51|). The scatter plot displays the 
values of £cr" c ^° n for the configurations of set M (part (a) of Fig. 15) and of set N (part 
(b)). It is most remarkable that a sizeable fraction of the configurations are accessible to 
investigation by direct detection, since the sensitivity in this kind of experiments is expected 



to be significantly improved in the near future [|S|. The dashed line in Fig. 15 shows the 
discovery potential in case of an improvement by a factor of 10 in current sensitivities, 
what is within reach in a short time. Our analysis shows an interesting interplay between 
experiments of direct search for particle dark matter and measurements of low-energy p's 
in space. This property would obviously be dramatically reinforced, should the indication 
about a possible annual modulation effect be confirmed by new data. In fact, it is very 
intriguing that many configurations of set M are indeed in the region singled out by the 
experiment of Ref. Q]. Finally, we notice that some configurations are actually excluded 
by the direct-search upper bound. This put emphasis on the potentiality of direct detection 
measurements in providing information on dark matter searches of different nature. 

Part (b) of Fig. 15 shows how the requirement of a sizeable contribution to the relic 
abundance makes somewhat thinner the set of configurations contributing to the highest 
values of £crf c "f^° n , but still leaves a significant number of configurations inside the closed 
region and, anyway, close to the current upper bound curve. Correlation between £cr™jj^ n 
and the neutralino relic density is given in Fig. 16. 
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VIII. SEARCH AT ACCELERATORS 



Let us turn now to the question of whether configurations of sets M and N might be 
explored at accelerators. 

LEP at y/s =192 GeV may explore the configurations with a neutralino mass up to ~ 50 



GeV |[52|| . Then, from Fig. 12 we see that LEP will be able to investigate only marginally 
the configurations of set M and N in the gaugino sector. Experimental investigation of 
larger masses requires future upgrading of the Tevatron or LHC For instance, TeV33 could, 
under favorable hypothesis, explore a range up to m x ~ 125 GeV [53]. In this case, all the 
higgsino configurations can be explored, as well as a large portion of the gaugino sector. 

A further illustrative point is offered by a scatter plot of set M in the plane m^-tan (3, 
displayed in Fig. 17 (m^ is the mass of the lightest CP-even scalar Higgs boson). The 
representative points of the set cover almost completely the Higgs physical region. Part of 
these supersymmetric configurations (the ones on the left side of the solid curves) will be 
explored by LEP at \/s = 192 GeV and at y/s = 200 GeV, with a luminosity of 200 pb _1 



per experiment [52] 



IX. CONCLUSIONS 

We have presented a new analysis of the cosmic-ray antiprotons flux, expected on the 
basis of secondary p's, generated by interactions of cosmic-ray primaries with the interstellar 
medium, and of a possible exotic primary source of p's, originated by neutralino-neutralino 
annihilations in the Galactic halo. 

Improvements over previous calculations of secondaries depend mostly on: i) the use of a 
two-zone propagation model for diffusion of cosmic rays in the halo instead of the standard 
leaky box model; ii) the inclusion of an energy-loss effect in the propagation properties of 
cosmic rays (important for the antiproton low energy range considered in this paper); iii) 
the use of the new data on primary cosmic-ray proton spectrum, as measured by IMAX [T^| 



and CAPRICE |20] 



The neutralino-induced p flux has been evaluated in a MSSM at the electroweak scale, 
which incorporates all current accelerator constraints. Use of supergravity-inspired unifica- 
tion conditions at large energy scale has been avoided in order not to arbitrarily constrain 



the neutralino phenomenology [53 



Solar modulation of the antiproton flux has been improved by analyzing the most com- 
plete set of data over the solar cycles p2[ and the data on the proton spectrum of Refs. 



HI- 

We have found that the most statistically relevant data on cosmic-ray antiprotons at low- 
energy |l[] leave some room for a possible signal from neutralino annihilation in the galactic 
halo. We have discussed how the relevant supersymmetric configurations may be explored 
with direct experiments for particle dark matter search and at accelerators. We have shown 
how the interplay between measurements of cosmic-ray p's and direct search experiments for 
relic particles is very intriguing and quite important in view of the significant improvements 
expected in these two classes of experiments in the near future. The present analysis stresses 
the great interest for the forthcoming AMS measurements with the Shuttle flight and on 
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the ISSA ||15|| , as well as for other future measurements with balloon-borne experiments 
(IMAX [|, BESS f|) and with satellites (PAMELA) gej, for disentangling the secondary 
p flux from a possible primary signal of exotic nature. As an example, we give in Fig. 18 
the distribution of measurements expected for AMS with the Shuttle flight according to two 
different hypothesis: a) dominance of the secondary contribution (lower sequence of crosses), 
b) significant contribution due a neutralino-induced signal (upper sequence of crosses). In 
our evaluation of the expected measurements we have taken into account geomagnetic cut-off 



effects and the expected AMS overall acceptance |55 
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TABLE I. Values of the parameters in the expressions (||) and (g) for the IS proton flux 
and of the solar-modulation parameter A. These values are obtained by best-fitting the data of 
Refs. [19-20] with Eqs. (|]) and (0), either over the entire energy range or only over the high-energy 
(T p > 20 GeV) range. First and third sets of values refer to 3-parameters fits (with Eqs.(|l]) and 
(|2|), respectively), second and fourth sets refer to 2-parameters fits at fixed A, (with Eqs.(|l]) and 
@, respectively). 





TA/T A V 

IIVIAA 




1 

Comments 


A 


12,300±1,700 


17,600±500 


entire energy 


a 


2.67±0.03 


2.81±0.01 


range 


A 


510±40 


390± 5 




A 


12,300±3,000 


19,600± 3,000 




a 


2.67±0.06 


2.85±0.04 


T p > 20 GeV 


A 


510 (fixed) 


390 (fixed) 




B 


16,200±2,000 


26,000±1,200 


entire energy 


7 


2.73±0.03 


2.91±0.02 


range 


A 


795±35 


710±10 




B 


13,700±4,100 


22,800± 3,700 




7 


2.69±0.06 


2.88±0.04 


T p > 20 GeV 


A 


795 (fixed) 


710 (fixed) 





Figure Captions 

Fig. 1 - TOA antiproton flux as a function of the antiproton kinetic energy. The 
experimental points are the BESS95 data jj]. The curves are the median (solid line), minimal 
(dotted line) and maximal (dashed line) secondary fluxes calculated in this paper, solar- 
modulated at the time of the BESS95 measurement. 

Fig. 2 - TOA spectra of IMAX (full circles) @ and of CAPRICE (open circles) @ 
with our best-fit curves with parametrization of Eq. ([I]). (The error bars are not shown 
when they are smaller than the dimension of the circles.) 

Fig. 3 - TOA spectra of IMAX (full circles) pf and of CAPRICE (open circles) |D 



The solid, dotted and dashed lines denote the median, minimal and maximal IS proton 
fluxes, respectively, as discussed in Sec. II. (The error bars are not shown when they are 
smaller than the dimension of the circles.) 

Fig. 4 - The grammage A e of the CNO primary elements (solid) as inferred from a 
two-zone diffusion model of the propagation of cosmic rays in the Galaxy. It is plotted 
as a function of the kinetic energy per nucleon. The dashed curve features the grammage 
corresponding to protons while the dotted lines delineate the interval of escape lengths 
inferred from the Ficenec et al. [|39| observations on 3 He at TOA energies comprised between 
100 MeV/n and 1.6 GeV/n. 

Fig. 5 - IS secondary antiproton spectra as functions of the p kinetic energy. Solid, 
dotted and dashed lines denote the fluxes obtained from the median, minimal and maximal 
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IS primary proton fluxes. The dot-dashed line denotes the median p spectrum, when the p 
energy losses are neglected. 

Fig. 6 - Coefficient C snsy (Tp, f) as a function of the p kinetic energy for different values 
of the flattening parameter /. 

Fig. 7 - Time variation of the solar-modulation parameter A. Full circles represent the 
best-fit values to the PGS average fluxes at minima (MIN) and at maxima (MAX) and to 
the fluxes of IMAX ^ and of CAPRICE J2g ; the open circle refers to the BESS95 data 
taking period and the cross denotes the extrapolated value of A at the time relevant for the 
future AMS Shuttle flight (May 1998). 

Fig. 8 - Solar modulation of the IS median secondary antiproton flux calculated in this 
paper. Solid line is the IS spectrum; dashed (dotted) line is the solar-modulated spectrum 
at minima (maxima). 

Fig. 9 - Solar modulation of the IS antiproton flux, due to neutralino annihilation for the 
representative neutralino configuration with m x = 62 GeV, P = 0.98 and f2 x /t 2 = 0.11. Solid 
line is the IS spectrum; dashed (dotted) line is the solar-modulated spectrum at minima 
(maxima). 

Fig. 10 - TOA antiproton fluxes versus the antiproton kinetic energy. The BESS95 data 
]T|] are shown by crosses. The dashed line denotes the median secondary flux, the dotted 
one denotes the primary flux due to neutralino annihilation in the halo for a neutralino 
configuration with m x = 62 GeV, P = 0.98 and Q x h 2 = 0.11. Solid line denotes the 
calculated total flux. 

Fig. 11 - Scatter plots for configurations of set M (a) and set N (b) in the P-tan/3 
plane. 

Fig. 12 - Scatter plots for configurations of set M (a) and set N (b) in the ~P-m x plane. 
Fig. 13 - Scatter plots for configurations of set R in the P-tan (3 plane (a) and in the 
P-m x plane (b). 

Fig. 14 - Scatter plots for configurations of set R in the P-tan j3 plane (a) and in the 
P-m x plane (b) for a flattening of / = 0.5. 

Fig. 15 - Scatter plot of the values of £a s n c u a f a e r on versus the neutralino mass for the 
configurations of set M (a) and of set N (b). The open curve denotes the (90 % C.L.) 
upper bound obtained from experimental data of Ref. p9 ]. The region delimited by a closed 
contour is the one singled out by the experiment of Ref. |50| as possibly indicative of an 
annual modulation effect. The total local dark matter density is normalized here to the 
value pi = 0.4 GeV cm" 3 . The dashed line shows the discovery potential in case of an 
improvement by a factor of 10 in current sensitivities for experiments of direct search for 
particle dark matter. 

Fig. 16 - Correlation between £cr" c ^ c / ar on and the neutralino relic density Q x h 2 for config- 
urations of set M. 

Fig. 17 - Scatter plot for configurations of set M in the m^-tan/3 plane. The region 
on the left of the dashed line denoted by (a) is excluded by current LEP experimental data 
P!f , the one on the right of the dashed line (b) is theoretically disallowed. The other lines 
display the LEP reach at luminosity L = 200 pb _1 and various energies |52|: (A) discovery 
potential at \fs = 192 GeV; (B) discovery potential at y/s = 200 GeV; (C) exclusion at 
y/s = 200 GeV. 

Fig. 18 - Expected distribution of measurements with the AMS Shuttle flight according 
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to two different hypothesis: a) dominance of the secondary contribution (lower sequence 
of crosses), b) significant contribution due a neutralino-induced signal (upper sequence of 
crosses). The dashed line denotes the secondary flux, the dotted one denotes the primary 
flux due to neutralino annihilation in the halo for a neutralino configuration with the repre- 
sentative values: m x = 62 GeV, P = 0.98 and Q, x h 2 = 0.11. Solid line denotes the calculated 
total flux. 
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